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Abstract 

We study the properties of fluctuation for the free energies and internal energies of 
two spin glass systems that differ for having some set of interactions flipped. We 
show that their difference has a variance that grows like the volume of the flipped 
region. Using a new interpolation method, which extends to the entire circle the 
standard interpolation technique, we show by integration by parts that the bound 
imply new overlap identities for the equilibrium state. As a side result the case of the 
non-interacting random field is analyzed and the triviality of its overlap distribution 
proved. 
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1 Introduction and results 



In this paper wc investigate a new method to obtain overlap identities for the spin glass 
models. The strategy we use is the exploitation of a bound on the fluctuations of a 
quantity that compares a system with some Gaussian disorder with the system at a flipped 
( J — > — J) disorder. While the disordered averages are symmetric by interaction flip due 
to the symmetry of the distribution, the difference among them is an interesting random 
variable whose variance can be shown to grow at most like the volume (for extensive 
quantities). 

The identities are then deduced using some form of integration by parts in the same 
perspective in which they appear from stochastic stability [AC] or the Ghirlanda Guerra 
method [GG] in the mean field case or, more recently, for short range finite dimensional 
models [CGi, CGi2] (see also [T, B]). 

The interest of obtaining further information from the method of the identities lies 
on the fact that they provide a constraint for the overlap moments (or their distribution) 
and have the potential to reduce its degrees of freedom toward, possibly, a mean field 
structure like it is expected to have the Sherrington Kirkpatrick model. 

More specifically the results of this paper consist of overlap identities for the quenched 
state which interpolate between a Gaussian spin glass and the system where the couplings 
in a subvolume (possibly coinciding with the whole volume) have been flipped. The inter- 
polation is obtained by extending to the whole circle the Guerra Toninelli interpolation 
[GT] . The bounds are derived from the concentration properties of the difference of the 
free energy per particle in the two settings, original and flipped. 

As an example, one may consider the result which is stated in [NS] (and quoted there 
as proved by Aizenman and Fisher) for the difference AF between the free energy of the 
Edwards- Anderson model on a d-dimensional lattice with linear size L and a volume L'^ 
when going from periodic to antiperiodic boundary conditions on the hyperplane which is 
orthogonal to (say) the x-direction. The mentioned property is a bound for the variance 
of this quantity which grows no more than the volume of the hyperplane. Such an upper 
bound is equivalent to a bound for the stiffness exponent 9 < {d — l)/2 [SY, BM, FH] (See 
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also the discussion of that exponent in [vE]). Although that bound is not expected to be 
saturated we prove here that it implies an identity for the equilibrium quantities. When 
expressed in terms of spin variables some of the overlap identities that we find generalize 
the structure of truncated correlation function that appear in [Te] whose behaviour in the 
volume is related to the low temperature phase properties of the model. Consequences of 
our bound can also be seen at the level of the difference of internal energies. This second 
set of identities contains as a particular case some of the Ghirlanda-Guerra identities. 

A quite interesting result, from the mathematical physics perspective, is provided by 
the analysis of the identities for the random field model without interaction. We show 
here that the new set of identities that we derive (and explicitly test) when considered 
together with the Ghirlanda Guerra ones provide a simple proof of triviality of the model 
i.e. the proof that the overlap is a non fluctuating quantity. We plan to apply the same 
method to the investigation of the random field model with ferromagnetic interactions. 

The plan of the paper is the following. In the next section we define the setting of 
Gaussian spin glasses that we consider. Then in section 3 we prove a lemma for the first 
two moments of the difference of free energies. This is obtained by studying a suitable 
interpolation on the circle for the linear combination of two independent Hamiltonians. 
Section 4 contains the proof of the concentration of measure results. The main results are 
given in section 5 and 6, where the new overlap identities are stated. Finally in section 7 
we study the case of the random field model and shows how to derive the triviality of the 
model without making use of the explicit solution. 

2 Definitions 

We consider a disordered model of Ising configurations (j„ = ±1, n e A C £ for some 
subset A (volume |A|) of a lattice £. We denote by Sa the set of all a = {cr„}„gA, and 
|Ea| = 21"^'. In the sequel the following definitions will be used. 

1. Hamiltonian. 

For every A C £ let {H\{a)}a-eSN be a family of 2l^l translation invariant (in 
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distribution) Gaussian random variables defined, in analogy with [RU], according 
to the general representation 

^A(a) = -J2 Jxax (2.1) 

XcA 

where 

(Tx=l[ai, (2.2) 
((70 = 0) and the J's are independent Gaussian variables with mean 

Av(Jx) = , (2.3) 

and variance 

Av(Jl) = . (2.4) 
Given any subset A' C A, we also write 

HA{a) = HAia) + HA\A,{a) (2.5) 

where 

HA'{a) = - ^ JxCTx, Ha\a'{(t) = - ^ Jxcrx , (2.6) 

XcA' ^CA 

X<t A' 

and 

^A,A'(^) = -HA'ia) + HA\Aia) (2.7) 

will denote the Hamiltonian with the J couplings inside the region A' that have 
been flipped. 

2. Average and Covariance matrix. 

The Hamiltonian Ha{(t) has covariance matrix 

CA(a,r) := Av {HA{a)HA{T)) 

= J2 ^\oxTx . (2.8) 

XcA 

By the Schwarz inequality 

XcA 

for all a and r. 



3. Thermodynamic Stability. 

The Hamiltonian (2.1) is thermodynamically stable if there exists a constant c such 
that 

sup ^ V < c < oo . (2.10) 

AC£ |A| 

Thanks to the relation (2.9) a thermodynamically stable model fulfills the bound 



CA{a,T) < c\A\ (2.11) 
and has an order 1 normalized covariance 

Ca((7,t) := i^CA((7,T). (2.12) 

4. Random partition function. 

Za(/3) := ^^"''''^''^ = Yl e-'^^A'('^)-'^^A\A'('^) , (2.13) 

2a,A'(/?) Yl e-^''^-^'^"^ = 5^ e'^^A'('^)-^^A\A'('^) . (2.14) 

ae Sa 0-e Ea 

5. Random free energy /pressure. 

-I3Tk{I3) := Vk{I3) := lnZA(/5) , (2.15) 

-/5^A,A'(/3) := n,A'(/3) := ln^A,A'(/3) • (2.16) 

6. Random internal energy. 

^a(/3) - ^^y- ' , (2.17) 

Z^o-e Sa 

E.^r i/AA'((T)e-^^A.A'M 

Z^CT6 Ea ^ 

7. Quenched free energy /pressure. 

-PFj,iP) := Pa(/3) := Av(7'a(/9)) . (2.19) 
-PFa,a'{P) := Pa,A'(/?) - Av(Pa,A'(/?)) ■ (2.20) 



8. R-product random Boltzmann-Gibbs state. 

[MP)] 



E (-) r^nmm ■ P-^D 



a(i),...,(T(«) 

9. Quenched equilibrium state. 

(-)a :=Av(Oa(-)) . (2.22) 

10. Observables. 

For any smooth bounded function G(ca) (without loss of generahty we consider 
IGI < 1 and no assumption of permutation invariance on G is made) of the covari- 

ance matrix entries we introduce the random (with respect to (— )) R x R matrix 
of elements {qk,i} (called generalized overlap) by the formula 

{G{q)) := Av(f7(G(cA))) . (2.23) 

E.g.: G(cA)=CA(a«,a(2))cA(a(2),a(3)) 

/ _ p-/3[E?=i^^a(^«)]\ 

(^1.2^2,3) = Av Yl ^a(^^'^ a(^))cA(a(^), a(^)) " (2.24) 

yiT(i),<T(2)^(j(3) y 

Remark: In the following, whenever there is no risk of confusion, the volume dependency 
in the quenched state or in the thermodynamic quantities will be dropped. 



3 Preliminary: interpolation on the circle 

Let ^ = {Ci}i<i<n a-nd rj — {r)i}i<i<n be two independent families of centered Gaussian 
random variables, each having covariance matrix C, i.e. 

Av(6e,-) = C,j 

MiiVj) = 0. (3.25) 
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Consider the following linear combination of ^ and rj 



where the parameter t e [a, 6] C R and the two functions f{t),g{t) take real values subject 
to the constraint 

ntr+g{tr = i. (3.26) 

Chosing f{t) — cos(i), g{t) = sin(t) we obtain: 

Xi{t) = cos(i) + sm{t) r]i. 

Because of the constraint (3.26), for any given time t e [a, b], the new centered Gaussian 
family x{t) — {xi{t)}i<i<n has the same covariance structure of ^ and rj: 

AY{Xi{t)Xj{t)) = dj, 

and hence the same distribution, independently of t (i.e. x{t) is a stationary Gaussian 
process) . 

In the abstract set-up described above, we regard x{t) as an interpolating Hamiltonian 
which is a linear combination of the random Hamiltonians C, and r], with t-dependent 
weights that are the coordinates of a point on the circle of unit radius.^ We introduce the 
interpolating random pressure ^: 

n 

V{t) = In Z{t) = In ^ e^'W ^ (3.27) 
1=1 

and the notation {€1^2) t,s to denote the expectation of the covariance matrix in the de- 
formed quenched state constructed from two independent copies with Boltzmann weights 
x{t), respectively x{s). Namely: 

„Xi{t)+Xj{s) 

(C,A« = Avi:C.,^^. (3.28) 



Ht is probably worth noting that any other parametrization of the unit circle would lead to the same 
expression as in (3.31). 

^Here, in defining the interpolating (random) pressure, we absorb the temperature in the Hamiltonian. 



The definition is extended in the obvious way to more than two copies. We will be 
interested in the random variable given by the difference of the pressures evaluated at the 
boundaries values 

X{a,b)^Vib)-Via) . (3.29) 

The following lemma gives an explicit expression for the first two moments of this random 
variable. 

Lemma 1 For the random variable X{a, b) defined above we have 

Av(A'(a,6)) =0 (3.30) 

and 

AY[{X{a,b)f] dtdsh{t,s){C,,2)t,s (3.31) 

J a J a 

- I I dtdsk2{t,s)[{Cl^)t,s-'^{C^,2C2,^)s,t,s+{Cl,2C^^^^^^^^ 
J a J a 

with 

ki{t,s) — cos{t — s), k2{t,s) — s\T?{t — s). (3.32) 

Proof. The stationarity of the Gaussian process x{t) implies that Av(P(t)) is independent 
of t, this proves (3.30). As far as the computation of the second moment is concerned, 
starting from 

Ay{X{a,b))^ j dtAy{r'{t))^ J dt Ay (x[it)^j (3.33) 

we have 

Ay[{X{a,b)f] = f dt f dsAy{V'{t)V'{s)) 

J a J a 

rb rb n , x,(t)^x,{s]\ 

= / */ E Av {^•m^')^^) ■ (3.34) 

The computation of the average in the rightmost term of the previous formula, which is 
reported in Appendix 1, gives 



Av \x^{t)x^{^s) 2{t)r^^s) ) ^ ~ «)(<^i,2)m 
sin^s - 1) (< C^,, -2 < C12C23 >t,s,t + < CuCu >t,s,s,t) (3.35) 



proving (3.31). 



□ 



4 Bound on the fluctuations of the free energy dif- 
ference 

It is a well established fact that the random free energy per particles of Gaussian spin 
glasses satisfies concentration inequalities, implying in particular self-averaging. Here 

we prove that the same result holds for the variation in the random free energy (or 
equivalently the random pressure) 

A-A.A' =Vk- Pa,a' (4.36) 

induced by the change of the signs of the interaction in the subset A' C A. In general, the 
fact that the random free energy per particle concentrates around its mean as the system 
volume increases of the free energy can be obtained either by martingales arguments 
[PS, CGi2] or by general Gaussian concentration of measure [T, GT2]. Here we follow 
the second approach. Our formulation applies to both mean field and finite dimensional 
models and, for instance, includes the non summable interactions in finite dimensions 
[KS] and the p-spin mean field model as well as the REM and GREM models. 
Before stating the result, it is useful to notice that, as a consequence of the symmetry of 
the Gaussian distribution, the variation of the random pressure has a zero average: 

Av(A'a,aO = . (4.37) 

Lemma 2 For every subset A' C A the disorder fluctuation of the free energy variation 
Aa,a' satisfies the following inequality: for all x > 

P(|^A,A.|>.) < 2exp(-5^) (4.38) 

with c the constant in the thermodynamic stability condition (cfr. Eq. (2.10)). The 
variance of the free energy variation satisfies the bound 

Var{X^^K>) = ^^{XIa) ^ 167rc/?2|A'| (4.39) 
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Proof. Consider an s > and let x > 0. By Markov inequality, one has 



P {A'a.a' >x} = F {exp[sA'A,A,] > exp(sx)} 

< Av(exp[sA'A,A']) exp(-sx) (4.40) 

To bound the generating function 

Av(exp[sA'A,A']) (4.41) 
one introduces, for a parameter t e [0, |], the following interpolating partition functions: 



Here the hamiltonians H^}{a), H^l){a), H^j^^^,{a), defined according to (2.6), depend 
on three independent copies {Jx^}xca, {Jx^}xca-> {Jx^}xca of the Gaussian disorder 
characterized by (2. 3), (2. 4). Now we are ready to consider the interpolating function 

<P{t) = ln^^;3A^;i jexp (^s Av^ |ln |^|^ | , (4.44) 

where Avi{—}, Av2{—} and Av^{—} denote expectation with respect to the three inde- 
pendent families of Gaussian variables Jx- It is immediate to verify that 

0(0) = In Av exp[s Aa,a'] , (4.45) 

and, using (4.37), 

</> (I) = . (4.46) 

This imphes that 

Av (exp[sA'A,A']) = e^^o)-'^^^) = e" <P'm_ (4.47) 



On the other hand, the function (f)'{t) can be easily bounded. Defining 

z-{t) 



K{t) = exp ( s Av2 { In |^ |. ) (4.48) 
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the derivative is given by 

+ (4-49) 

where 

^ ^^_ sAv,Av,{K{t)Av,{^}} sAvsAv,{K{t)Av,{Pi:^}} 
Av^Av,{K{t)} ' AvsAv,{K{t)} 

(4.50) 

The first term in the derivative is 

sAv^Av^ {K{t)Av2 {j:,e^^pti<^) [/3 sin i //jV (a) - (3 cost H^^) (a)]}} 

where 

-/3 cos t if^V (t) -/S^' A' * ^A'^ 

"'^^^ = ZHT) 

Applying the integration by parts formula, a simple computation gives 
psintAvsAvi Av2 |^p+((T)i/iV(a) 

= -s/^^sint costAv^Avi I K{t) ^ A|[s+(X)2 + s+(X)s7(X)] 

I XcA' 

- /S^ sin i cos i Av^Avi | | J]] Ca/ (cr, (7)pt+ (a) 1 1 

+ P'^sint costAv^Av^ \^K{t) Av^ |^Ca'((7,t)p+((7)p+(t))| | (4.53) 



and 



-PcostAv^Av^ yHt) Av2 <^p\(t) ^^l^V) 
= sin t cos t Av^Avi | K{t) Av2 | ^ Ca' (o", o-)Pt^ } } 
- P^sint costAvi I K{t) Av2 \ Ca^^^, r)p+((7)p+(r)) \ \ (4.54) 



cr,T 



where 



st{X) = Av2 |x]axp+(a)| , s-{X) = Av2 ^axP^{a)^ (4.55) 
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and 

J cos t H';> {T)-PH)^^^, sint H)^/ {a, 

= z^) ■ (^•^«) 

Taking the difference between (4.53) and (4.54) one finds that 

m ^ -.y sint cost ^Xt'gwf ^ ^' . (4.57) 

With a similar computation one obtains also 

<t>_{t) = /5 smt cost Av,Av,{m] ' ^'-''^ 

then we conclude that 

m - -s^f^^t COS* ^"''^"^ ^:x;^Sr' ^ ■ 

Using the thermodynamic stability condition (2.11), this yields 

\(t)\t)\<A(5^cs^\h'\ (4.60) 

from which it follows, using (4.47) 

Av (exp[sA'A,A']) < exp {2t:(5'^cs^ |A'|) . (4.61) 
Inserting this bound into the inequality (4.40) and optimizing over s one finally obtains 

P(Af,,.>.) <exp(-^^^). (4.62) 

The proof of inequality (4.38) is completed by observing that one can repeat a similar 
computation for P (A'a^a/ < —x). The result for the variance (4.39) is then immediately 
proved, thanks to (4.37), using the identity 

Av =2 X P(|Aa,a'| > x) dx . (4.63) 

Jo 

□ 
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5 Overlap identities from the difference of free en- 
ergy 

We are now ready to state our first result. 

Theorem 1 Given a volume h., consider the Guassian spin glass with Hamiltonian (2.1). 
For a suhvolume A' C A and a parameter t e [0, tt], let 

a 

with 

H,{t) = cos{t)H^l){a) + sin{t)Hf){a) + H^\^,{a) 

he the Boltzmann-Gibbs state which interpolates between the system with Gaussian dis- 
order and the system with a flipped disorder in the region A' ( Hj^} and Hf) are two 
independent copies of the Hamiltonian in the subvolume A', H\\\>{a) is the Hamilto- 
nian in the remaining part of the volume, they are all independent). Then, the following 
identities hold 

lim / dtds sm\s - t) {{c^',)\s - 2(<2<3).,m + {ci,2C3A)t,s,s,t = (5.64) 

where ((cf'g)^)*,^ (and analogously for the other terms) is the overlap of region A' Q A in 
the quenched state constructed form the interpolating Boltzmann-Gibbs state, e.g. 

((<2)')t,. = Av(a;ta;,(4,(a,T))) . 

Proof: The proof is obtained from a suitable combination of the results in the previous 

sections. For a parameter t G [0,7r] we consider the interpolating random pressure 

V{t) = In ^ e^-(*)+-^MA'W (5.65) 

where 

x^it) = cos(t)i/J'(a) + sin(t)i/|'V) 

with H^}{a), H^}{a) two independent copies of the Hamiltonian for the subvolume A' C 
A. The boundaries values give the random pressure of the original system when t — 
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and the random pressure of the system with the couphngs J flipped on the subvolume A' 
when t — TT, i.e. 

V{0) = Va , 

Apphcation of Lemma 1 with = H^} (a) and rja — H^} (a) (the presence of the addi- 
tional term HA\A'i<7) in the random interpolating pressure does not change the result in 
the Lemma, as far as the quenched state is correctly interpreted) gives 

/•TT /"TT 

VariVA-VA,A') = A' / / dt ds cos{s - t){4'2)t,s (5.66) 

Jo Jo 

+ (A')' £ [ dt ds sin^(. - t) [{{4;,)\s - 2{ct;,4;,)s,t,s + {ct24'As,s,t] 

On the other hand, Lemma (2) tell us that Var{VA—VA,A') is bounded above by a constant 
times the subvolume A'. As a consequence, the statement of the theorem follows. □ 
Remark: When expressed in terms of the spin variables the polynomial in the integral 
(5.64) involves generalized truncated correlation functions. Indeed, for the model defined 
in Section 2, we have the following expressions 



X,YCA' 

^s,tM',2C^',^ = E A|A^a;,(a«)a;,(aJVgVs(4'^) 

X,YcA' 

u;t,s,sAC^:2C^:d = Yl A|A^u;,(a«)a;,(aJVs(4'Vt(4'^) (5-67) 

X,YcA' 



thus 



^m((<2)') - 2 Ws,m(<2C2'3) + ^M,s,t(<2<4) = (5-68) 

Y A^Ay [utiaxcry) - <^t((7x)a;t((7y)] [a;s((7x(7y) - u;s{crx)uJs{crY)] 

' ' X,YCA' 

where replica indices have been dropped. For the Edwards-Anderson model, which is 
obtained with A^ = 1 if X e = {{n, n') e A' x A', \n-n'\ = 1} and A^ = otherwise, 
the linear combination (5.68) of the moments of the link-overlap in the region A' is written 
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in terms of truncated correlation functions, that is 

^ [ut{abai,^) -uJt{ab)uJt{ab,)][u}s{(^b(^b') -i^s{crb)i^s{crb')], (5.69) 

' ' b,b'eB' 

with (jft = CnC^, a b = {n, n') e B' . 

6 Overlap identities from the difference of internal 
energy 

In this section we study the change in the internal energy after a flip of the couplings. 

We consider only the case of the flip of all the couplings in the entire volume. 

Let us consider two centered gaussian families ^ — {^i}i<i<n, ^ = {Vi}i<i<n with 
covariance structure given by 

Av{Uj)=Av{7^i7^j)=Ci,j (6.70) 

with Cj.i = A^. We assume the thermodynamic stability condition to hold. It follows 
that N is proportional to the volume. For example, in the case of the Edwards-Anderson 
model on a o?-dimensional lattice we would have N — d\A\. We introduce the random free 
energies 

Vd(3)^ In Z^{(3)^ In P,(/?) = lnZ,(/?) = In^e'^^S (6.71) 

i i 

with the random Boltzmann-Gibbs state a;^(— ), a;^(— ) and their quenched versions: 

(-)^ = Av4(^^(-), (-), = Av,(^,(-). (6.72) 

With a slight abuse of notation we will use the previous symbols also to denote the product 
state acting on the replicated system. The free energy difference, obtained flipping the 
hamiltonian 77, 

= - ^ In ^ e-^«^ - In ^ e^^% (6.73) 

i i 
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has a /3-derivative given by the difference between the internal energies: 



Using the symmetry of the distribution of 77, we have the identities^ 

Av^a;^(0 = -(3{N - Av^a;^(Ci,2)) = -(3{N - (Ci,2)^) 
Av,a;_,(r;) = /3(iV - Av,a;_,(Ci,2)) = ^{N - (^1,2),). 

The above formulae show that the disorder average of X'{I3) vanishes 



(6.74) 



(6.75) 
(6.76) 



Av^,,(A"(/5)) = /5((Ci,2), - (Ci,2)«) = 0, 



(6.77) 



since, obviously, (Ci,2)?j = ((^1,2)^- Here Ci,2 = represents the covariance matrix 

whose entries are regarded as configurations of two replicas labeled 1 and 2. Thus, using 
the identity Av^(a;_^(r;)^) = Av^(a;^({)^), we have that the variance of A"(/5) is given by: 



ky^,,{X\(3Y) = 2ky^{uj^{iY) + 2Av^,,(a;^(0c^-,(r;)). 



(6.78) 



Using the integration by parts formula, we obtain that 



Av^(u;^(0^)=Av^J]a,- 
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. (6.79) 



The second term in the right-hand side of the previous formula requires a repeated appli- 
cation of the integration by parts formula, which gives: 



i,3 k,e 



-P'N{N-2{C,,,)^)+(3\Ci,)^ 

-6/3'(Ci,2C2,3)^ + 6/3^(Ci,2C3,4)e 



^Indeed, from the symmetry of the gaussian distribution, we have that for any function f{r]) the 
following equalities hold: Av^ifiv) = ^'^rif{—v) = Av_,,/(— 77). In particular if 5 is a function of the 
configurations of the replicated system, applying the previous remark to f{r]) = w^(fif) we obtain: (g)^ = 
AVrfOJj^{g) = AVrjUJ^rjig) = Av_^Lj_^(g) = {g)-r,. 

These properties will be tacitly used several time in this section. 
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Since the first term in the right-side of (6.79) is quenched average of Ci,2, we conclude 
that 



Ayd^dO') = {C,,2)+P'N{N-2{C,,,)) + p'{Cl) 

- 6/32(Ci,2C2,3> + 6/3^(Ci,2C3,4) (6.80) 

dropping, here and in what follows, the unessential reference to ^ in the quenched averages. 
If the two families ^ and 77 were independent, then in (6.78) the average of the product 
would factorize Av^,^(a;^(^)a;_^(77)) = —P'^{N — (Ci,2))^ giving: 

Av^,,(A"(/3)2) = 2(Ci,2) + 2(3' ((Cj^) - (Ci,2)') + 12(3' {{C^^^Cs,^) - {C,,2C2,s)) ■ 

In this case the self averaging of the normalized quantity X'{(3)' /N (see Theorem 2) would 
lead, in the large volume limit A?^ — > 00, to the well known identity [G] 

(C1,2C2,3) - (C1,2C3,4) = \ {{cl^) - {c^,2?) ■ (6.81) 

However, our concern here is the computation of the quadratic fluctuations of X'[(3) when 
the sign of a given hamiltonian ^ is flipped in the whole volume. Therefore we have to set 
^=77 in (6.78). The computation requires, once again, the repeated use of the integration 
by parts formula 



& 



2 



Z^{(3)Z^{-(3)^ 



(6.82) 



The average in (6.82) is expressed through a set of mixed quenched state: for instance, 
the first term in right-hand side of the previous equation is 

(C,,),,- = Av,^C„^j^^^^. (6.83) 

Generalizing the previous definition we have, for instance, that (— )+^+^_^+ represents 
the thermal average taken with the usual boltzmannfaktor (i.e. with the sign — in the 
exponent) in the first, second and fourth copy, and with the opposite sign in the third 
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one. Moreover, the symbol (— with all the subscripts + (or — , because of the 
symmetry of the gaussian distribution), is the usual quenched measure ( — ). The explicit 
computation gives: 

92 



Z^{(3)Z^{-(3) 



-4/3^(Ci,2C2,3)+,+,- + 4/32(Ci,2C3,4)+,+,+,- - /3^(Ci,2C3,4)+,+,-,- - /3^(Ci,2C3,4)+,-,+,- 

and finally: 

Av,(A"(/3)2) = 2((Ci,2)+,+ + + {{Clh,+ - (6.84) 

- 4/3^ (3(Ci,2C2,3)+,+,+ - (Ci,2C2,3)+,-,+ + 2(Ci,2C2,3)+,+,-) 

+ 2/3^ (6(Ci,2C3,4)+,+,+,+ + 4(Ci,2C3,4)+,+,+,- - (Cl,2C3,4)+,+,-,- " (Cl,2C3,4) +,-,+,-) • 

If we choose now ^ to be the Hamiltonian family defined in section 2, we obtain the 
following: 

Theorem 2 Consider the Guassian spin glass with Hamiltonian ^ given in (2.1). In the 
infinite volume limit and for almost all values of (5, we have 

[(c?,2)+,+ - (C?,2)+ -] - 2 [3(ci,2C2,3)+,+,+ - (ci,2C2,3) +,-,+ + 2(ci,2C2,3)+,+,-] 
+ [6(ci,2C3,4)+,+,+,+ + 4(ci,2C3,4)+,+,+,- - (01,203,4)+,+,-,- - (01,203,4)+,-,+,-] = (6.85) 

where (01,2)+,- (and analogously for the other terms) is the overlap expectation in the 
quenched state constructed form the mixed Boltzmann-Gibbs state with one copy given by 
the original system and the other copy given by the flipped systems, e.g. 

{4,2)+,- = Av(a;^a;-^(o^((7,T))) . 

Proof. 

The proof is a simple consequence of well known results. The sequence of convex functions 
V(^{l3)/N converges almost everywhere in J to the limiting value a{[3) of its average and the 
convergence is self averaging (i.e. Vax{V^{j3)/N) 0). By general convexity arguments 
[RU] it follows that the sequence of derivatives V'^{I3)/N converges to u{l3) — a' {(3) al- 
most everywhere in (3 and also that the convergence is self averaging (Va,-r{V'^{P)/N) — > 0, 
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/3-a.e.) [S, OTW]. These remarks apply obviously also to V-^{P)/N and to its deriva- 
tive, with the same limiting functions a(/9) and a'(/?). Thus we have that X[(3)/N — 
V^{P) /N—V-^{P) /N and its derivative X'{P) /N vanish a.e. in J in the large volume hmit. 
Moreover, Var(A"(/3) /N) = Var(7'^(/3) /N)+Ysij:(V'_^((3) /N) -2cov {r'^{(3) /N, VL^{/3) /N) , 
thus estimating the covariance with the Cauchy-Schwartz inequality we have 

Var(A"(/?)/Ar) < 4Var(:P^(/?)/Ar) ^ 0, (3 - a.e. (6.86) 

for — > oo. Therefore, dividing (6.84) by N'^ and taking the limit we obtain the result. 

□ 

7 Triviality of the Random Field model 

In this section we compute explicitly the expression appearing in Theorem 1 

(Cl,2)t,s - 2(ci,2C2,3)s,t,s + {Ci^2C^,i)t,s,s,t (7.87) 

in the simple case of the random field. We will show that this linear combination of 
overlap moments vanishes pointwise for all values of t and s. We will then deduce the 
triviality of the order parameter for the random field model. 

We consider two families Jj and Jj for i = 1, A?" of independent normally distributed 
centered random variables with variance 1: 

Av(JiJ,) = Av(JiJ,) = Av(JiJ,) = 0, (7.88) 

and the random field hamiltonians 

N N 

^a = ^ JiCTi, 'na = ^ Jm- (7.89) 

i=l i=l 

where CTj — ±1. We have that ^ = {^aja and r] = {ri^j}^ are two independent centered 
Gaussian families (each having n = 2^ elements indexed by the configurations cr, N being 
the volume) and covariance structure given by: 

A.v{l]aVT) = Ca,T = Nq{a,T), 

MCaVr) = 0. (7.90) 
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where q{(J,T) is the site overlap of the two configurations a and r: 

1 ^ 

The interpolating Hamiltonian: 

x^{t) = cos(t)C^ + sin(t)77^, (7.92) 
which is a stationary Gaussian process with the same distribution of ^ and rj: 

Av{x,{t)xr{t)) = Nq{a,T), (7.93) 

defines the quenched deformed state on the rephcated system, whose averages are denoted 
with the usual notation, e.g. {—)t,s, {—)s,t,s ■ ■ ■ ■ 

Theorem 3 Consider the random field spin glass with Hamiltonian (7.89). In the limit 
— > oo and for all values of t and s we have 

ll{(ll,2)t,s + l2{qi,2)t,s + I3{qi,2q2,3)s,t,s + M(li,2q3,4)t,s,s,t = (7.94) 

for any choice of real 71, 72, 73, 74 with 71 + 72 + 73 + 74 = 0. 

Proof The simple proof relies on the following identities, derived in Appendix 2: 



{Cl,2)l 


N 

= ^(Av( tanh(G',(i)) tanh(G',(s))))2 + Qjv(i, s), 

i=l 


(7.95) 


1,2) t,s 


= l + QN{t,s), 


(7.96) 


(Ci,2C2,3)s,t,s 


N 

= J]Av(tanh2(G',(s))) + Q^(i,s), 
1=1 


(7.97) 


{Ci2C^4)t,s,s,t 


N 

= ^Av (tanh2(G',(t))tanh2(G^.(s))) + Q^(t,s), 


(7.98) 



where 

Gi{t) = cos{t)Ji + sm{t)Ji, (7.99) 
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and QN{t, s) is a term of order N^, see (9.118). Thus: 

7l(Ci%)t,s + 72(Cl,2)t + l3{Cl,2C2,3)s,t,s + -f3{Cl,2Cs,i) t,s,s,t = 
N N 

71 + 72 ^(Av( Unh{Gi{t)) Unh{Gi{s))))^ + 73 J] (tanh2(G'i(s))) 
1=1 1=1 

+74 J] Av (tanh2(G'j(i)) tanh^(G'j(s))) 

j=i 

+ (71 + 72 + 73 + l4)QN{t, S), 

i.e. the hnear combination of the covariance matrix moments is of order A^. Thus, since 
I tanh(a;)| < 1, we have 

|7i(Ci,2)m + 72(Cl,2)J, + 73(Cl,2C2,3)s,M + l4{Cl,2C3,i) t,s,s,t\ 

< |7i| + (I72I + I73I + |74|)A^, 
which can be rewritten, using the overlaps gi,2, ^2,3, ?3,4 between rephcas, as 

|7i(9?,2)m + 72(9l,2)t% + 73(gi,2g2,3).,M + 74 (gi,2g3,4)M,^,t | 

^ I72I + I73I + I74I |7l| rr. .^^^ 

□ 

Among the relations of theorem 3, in the thermodynamic limit, we find the identity of 
theorem 1 for the values 71 = 1, 72 = 0, 73 = —2, 74 = 1: 

(?l,2)t,s - 2(gi,2g2,3)s,t,s + {qi,2q3,4)t,s,s,t = (7.101) 

and the Ghirlanda-Guerra identities: for 71 = 1, 72 = 1, 73 = —2, 74 = we find 

for 7i = 1, 72 = 2, 73 = 0, 74 = —3 we find 

1 2 

(91,293,4).,*,. = 3(91,2)*,* + 3(91,2)?,. • (7.103) 
Using (7.102) and (7.103) we can express (7.101) as: 

(91,2)*,. - 2(91,292,3).,*,. + (91,293,4)*,.,.,* = ^((91,2)*,. - (9i,2)t,J (7.104) 
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The identity derived from the flip of the couphng thus imply a trivial order parameter 
distribution. Indeed, since the identity (7.101) is true for every t and s we can choose 
t — s — and then the interpolating states reduce to the usual quenched Boltzmann-Gibbs 
state. Prom Eq. (7.104) we deduce a trivial overlap distribution. 
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8 Appendix 1 

In this appendix we will use the Gaussian integration by parts formula for correlated 
Gaussian random variables zi, . . . ,Zn: 

Ay{zM^u . . . , ;^n)) = Av {z,z,) Av ^ (8.105) 

to compute the second moment of the pressure difference X{a,b). We have to evaluate 
the average inside the integral (3.34) 

n 

J2 Av {x[{t)xr{s)B{t,j;t,s)) = (8.106) 

n n 

= sm{t) sin(s) ^ Av {^i^jB{i,j; t, s)) - sm{t) cos(s) ^ Av {^ir)jB{i,j; t, s)) 

i,j=l i,j=l 



n 



- sin(s) cos(i) ^ Av {^jriiB{iJ; t, s)) + cos(i) cos(s) ^ Av {riirijB{iJ; t, s)) 

i,j=l i,j=l 



22 



where, for the sake of notation, we have introduced the symbol 



B(i,j;t,s) = 



Z{t)Z{s) • 
Applying (8.105) twice, we obtain 

Ay{^i^jB{i,j;t,s))^CijAy{B{i,j;t,s))+J2 C^,kCj,tA^ i-^B{i, j;t, s) j (8.107) 

k,e=i ^ ^'^^^ ' 

Av(7?i77,S(i,j;i,s)) =CijAv(S(i,j;i,s))+ C^.^C.-^Av (^^_-S(i, j; s) j (8.108) 

/ ^2 \ 

Av(ei?7.S(i,j;i,s)) = Av(e,-77iS(i,j;i,s)) = ^ Ci,feC,-,Av (^^S(i, j; s)j (8.109) 

The combination of the first two terms in the right hand sides of (8.107) and (8.108) 
with the trigonometric coefficients given by (8.106) produce the quenched expectation 

COs(i - s)(Cl,2)i,s- 

The explicit computation of the derivatives is long but not difficult; the result is: 

^^^^B{i,j; t, s) = B{i,j; t, s) {cos'^{t)Ai + cos'^{s)A2 + cos(t) cos(s)(A3 + A^)} , 
^^^B{i,j- 1, s) = B{i,j; t, s) {sin2(s)Ai + sin2(s)A2 + sin(i) sin(s)(A3 + ^4)} , 
^^B{i, j; t, s) = B{i, j; t, s) {sin(t) cos(t) Ai + sin(s) cos(s) 
+ sin(t) cos(s)743 + sin(s) cos(t)744} , 

where Ai, A2, A^, A4 are combinations of Kronecker delta functions depending on the 
indices i,j,£, k and Boltzmann weights for the hamiltonians x{t) and x{s). 
Using the previous formulas for the second derivatives and formulas (8. 107), (8. 108) and 
(8.109), we see that the right hand side of (8.106) contains a hnear combination of func- 
tions Aj with trigonometric coefficients given by the product of four factors taken from 
{cos(t), sin(t), cos(s), sin(s)}. It is not difficult to recognize that the coefficient of A^ is 
— sin2(s — t) while the other are zero. Thus, 



n n 



J2 Av {x[{t)x'j{s)B{i,j;t,s)) = cos(t - s){Ci^2)t,s - sm\s - t)Av Ci,kCj,eAsB{t,j;t,s) 

i,j=l i,J=l k,e=l 
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and since 

As = 5e,A, - - Skj^ + -z^Yis) 

we obtain 

n n 

Av ^ ^ Ci,kCj,iA3B{i,j;t,s) =< >t,s -2 < C12C23 >t,s,t + < C12C34 
i,j=i k,e=i 

which proves (3.35). 



9 Appendix 2 

In this appendix we prove the identities (7. 95), (7.96), (7. 97), (7.98). Recalhng the definition 
(7.99) of the Gaussian variables Gi{t), we can define the interpolating partition function 

Z{t) = E.exp(a:.(t)) = E. exp(EL G.(t)a.). (9.110) 

A simple computation shows that 

N 

Z{t) = 2^]^coshGi(t). (9.111) 

i=l 

For any integer M < N, we consider the sublattice Am = {N — M + 1, . . . , N} C 
Atv = {1,...,N} with its spin-configuration space Sm = { — 1,1}''^*'' and the subspace 
— {(+1) ■ ■ ■ , Cat), CTi — ±1} (we will drop the subscript when it is equal to 

N, e.g. S = Sn, = S^, etc.). The interpolating Boltzmann-Gibbs random state on 
the lattice A^ is 

^ Eae5M/(^)^^pfc=iV-M+lG'i(t)(Xi) 
= ^) ' 

ZM{t) = ^ cxp ( Yl G,{t)aA , (9.112) 

aeSM \i=N-M+l I 

where / is a function on Su- This definition extends in the obvious way to the i?-fold 
product; for instance the 2- product measure for the parameter values t and s is given by 

Ea,re5M /('^' ( E^=JV-M+l Gi{t)Ui + E^IjV-M+l Gi{s)T^ 

^tiif) = ■ ^ ^ , , , ^ 9.113 

'''^■^^ ZM{t)ZM{s) ^ ' 
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where / is a function on Sm x Sm- In the sequel we will also write ujt,s{—) instead of 

The computation of the moments of the covariance matrix is done evaluating (by induction 
on M) the averages of the products of the overlaps between configurations of Sm 

1 ^ 

i=N-M+l 

Indeed, the explicit computation shows that: 

= ^tanhGi(i)tanhGi(s) + ^L^u;^-\q^_,), (9.115) 
then iterating the previous formula — 1 times we obtain 



1 ^ 

<M = ^ $]tanh(G,(t)) tanh(G,(5)). (9.116) 



since ul ^{qi) = tanh(G'Ar(t)) tanh(Giv(s))- Recalling the relation between overlaps and 
covariance (7.90) and taking the average with respect to the disorder, we obtain: 

N 

{Ci,2)t,s = ^Av(tanh(G',(t))tanh(G'i(s))), (9.117) 
1=1 

thus 

N 

{Ci,2)l = ^(Av(tanh(G',(i)) tanh(G',(s))))^ + QN{t,s) 



i=l 

where 

Q^it, s) = 2 Av [tanh(Gj(t)) tanh(Gj(s))] Av [tanh(G^(t)) tanh(G^(s))] (9.118) 

i<j<e<N 

is a term of order A^^. This proves (7.95). 

For the squared overlap the following relation holds 

1 2 

^tA<lN) = ]7 + ]^ - tanh(G,(t)) tanh(G,(.))a;fr^(g^_i). (9.119) 



Since for M < iV 

N 

M 



1 ^ 

<s(9m) = t7 E tanh(G,(i))tanh(G,(5)) (9.120) 



j=N-M+l 
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we can write 



^tMl) = T72+^ tanh(G,(t)) tanh(G,(s)) tanh(G,(t)) tanh(G,(s)) (9.121) 



7V2 Ar2 

i<j<^<Ar 



and finally 



{Cl^)t^, = 1 + 2 Av [tanh(Gj(t)) tanh(Gj(s)) tanh(G^(t)) tanh(Gf(s))] . 

1<3<1<N 

From the independence of the random variables Gi{t) (see (7.99)), we have that the 
average in the right hand side of the previous formula factorizes, thus we obtain (7.96): 

(C?,2)t,. = l + Qiv(t,s)- 

The second term in (7.87) is computed considering the average of ?iv(cr, 7)g7v(7, t) where 
7, (7, T e S. We have 

'^i!t,.(?iv((T,7)giv(7,r)) = -^tanh2(G'i(s)) 
+-^tanh(G'i(i)) tanh(G'i(s)) Y!^=2'^^^HG jit)) tanh(G'^(s)) 

+ (^)'<,-^(giv-i(a',7')?iv-i(y,r')). (9.122) 

where cr', 7', r' are the restriction of cr, 7, r to Sn-i- As in the previous cases, iterating this 
formula and taking into account that a;^j^^(gi((T, 7)^1(7, r)) = tanh^(Gjv(-s)) we obtain 

'^fM(9^(^,7)9iv(7,r)) = ^Eiitanh^l^ils)) 
+^ Ef="i' tanh(G', (i)) tanh(G',(s)) Ef=,+i tanh(G',(t)) tanh(G',(s)), (9.123) 

then 

TV 

(C^i,2C2,3).,M = (tanh^(G',(s))) + QN{t,s). 



i=l 



which proves (7.97). The computation of the last term in (7.87) is simple because in this 
case the random product measure factorizes: 

(^t,s,s,Mi.^^ '^)5(7, «^)) = ^m(?(o-, T))uJs,t{q{l, «;))■ (9.124) 
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Then, using (9.116) we have 

1 ^ 

wt,,,,,t(g(a,r)g(7,«:)) ^ j^Yl tanh(G',(t)) tanh(G',(s)) tanh(Q(t)) tanh(Q(s))(9.125) 
and 



N 

{Ci2C34)t,s,s,t = 5^ Av (tanh(G',(t)) tanh(Gj(s)) tanh(G'i(t)) tanh(Gi(s))) (9.126) 

which, using the symmetry of Ojj = Av (tanh(Gj(t)) tanh(Gj(s)) tanh(Gj(t)) tanh(Gj(s))), 
gives (7.98): 

N 

(Ci2C34)t,.,.,t = 5]Av (tanh2(G,(t))tanh2(G',(s))) + Qiv(t,s). 
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